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Abstract 

We have analyzed the tunneling transmission probability and electronic cur- 
rent density through resonant heterostructures in the presence of an external 
electromagnetic field. In this work, we compare two different models for a 
double barrier : In the first case the effect of the external field is taken into 
account by spatially dependent AC voltages and in the second one the elec- 
tromagnetic field is described in terms of a photon field that irradiates homo- 
geneously the whole sample. While in the first description the tunneling takes 
place mainly through photo sidebands in the case of homogeneous illumination 
the main effective tunneling channels correspond to the coupling between dif- 
ferent electronic states due to photon absorption and emission. The difference 
of tunneling mechanisms between these configurations is strongly reflected in 
the transmission and current density which present very different features in 
both cases. 

In order to analyze these effects we have obtained, within the Transfer Hamil- 
tonian framework, a general expression for the transition probability for co- 
herent resonant tunneling in terms of the Green's function of the system. 
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I. INTRODUCTION 



In the last years several works have been devoted to the analysis of the effect of a time- 
dependent field on the transport properties of resonant heterostructures, i.e., double barriers, 
quantum wires and quantum dotffl. However there is not yet a systematic discussion of 
the different situations corresponding to an AC voltage applied between the left and right 
leads and which implies a position-dependent dephasing of the external field on the sam- 
ple , and the case where the whole sample is homogeneously illuminated ,i.e. , where the 
electron-photon coupling depends on position just through the momentum matrix elements. 
In spite of the increasing interest in this field, most of the theoretical work has been done 
considering spatially time dependent voltages but the experimental information includes 
both configurations. Both situations, however, are physically different because the effective 
tunneling channels for the flowing of current are different in both cases. 
The coherent tunneling in the presence of light applied homogeneously on a double barrier 
heterostructure has been recently treated in the scheme of the Transfer MatrixS. 
An alternative way to obtain the Transmission coefficient is based in the Transfer Hamilto- 
nian and considers a localized basis representation, i.e., approximated hamiltonians whose 
eigenstates are spatially localized. The extension of the TH to analyze coherent resonant 
tunneling (GTH)E-a allows to analyze not only the sequential tunneling which consider the 
electrons tunneling through the single barriers, emitter and collector, sequentially, but the 
coherent one which includes the virtual transitions through the resonant states for electrons 
crossing coherently the heterostructure. 

In this work we have extended the Generalized Transfer Hamiltonian formalism (GTH)0 
to obtain the transition probability for the coherent tunneling in the presence of a time 
dependent potential. 

In the first section of the paper we will develop the theory to obtain the transition proba- 
bility for a double barrier structure in the presence of a spatially dependent AC modulation 
. In the next section the GTH is extended to describe the photoassisted tunneling process 
for a sample homogeneously illuminated . The different features obtained for the coherent 
tunneling current for the two different configurations are discussed in a further section. 



II. TRANSITION PROBABILITY THROUGH A DOUBLE BARRIER IN AN AC 

FIELD 

We are going to analyze the effect of an AC field applied just to the left and right leads 
with a dephasing of n . This configuration is schematically represented for the double barrier 

Th^TVansfer Ha—, fetalis, developed by BardeeaB , allows to describe in & st 
order time dependent perturbation theory the transition probability in terms of the eigen- 
states of auxiliar hamiltonians spatially localized. This formalism has been extended to all 
orders in perturbation theory (Generalized Transfer Hamiltonian) to analyze as well tunnel- 
ing through resonant states, i.e., to include virtual transitions through the localized states 
in the well" (see fig. 2a). 
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The quantum mechanical Hamiltonian for an electron in the presence of time dependent 
potential can be written: 



Hit) = H L (t) + H R (t) + H T (t) + H.c (1) 

Where Hi(t) and H R (t) are the hamiltonians for the left and right sides respectively includ- 
ing the time dependent perturbation and Hx{t) is the coupling term which accounts for the 
transfer of electrons from the left to the right side. We are going to define now H^t) and 
H R (t) as: 

H L (t) = X>z.(*)<c fci = J2[e kL + (k L \H 2 (t)\k L )]4 L c kL 

k k 

H R(t) = E e PflW c p« c P« = J21 € pr - (PR\ H 2(t)\PB)]Cp R c PR . (2) 
p p 

where c£ L , Ck L , c+ fl , c PR are creation and destruction operators for electron in the left- and 
right-hand side of the heterostructure, respectively; efc L (e Pfl ) are the eigenenergies for 
(H R ) and H 2 (t) is the hamiltonian which describes the AC time dependent potential in the 
leads and can be written as: 

H 2 {t) = V Ac cosu t (3) 

In these expressions just diagonal terms are considered in the time dependent term. The 
reason is that the off-diagonal terms are zero in this case due to the orthogonality of the 
eigenstates of the auxiliar hamiltonians H®,H R . The retarded single electron Green's func- 
tion of these hamiltonians are: 



G+(M') = -^0(t-t')(0\c kL (t)4 L (f)\0) = - l -e(t-t')exp[- % -J ti d^e.M] 

■ oo oo \r 1/ 

G + R {t,t') = -^e(t-t'){0\c PR (t)c; R (t')\0) = -^(t-t^expi-'-j^dhe^h)] 

• OO OO T/ \T 

= ~i E E M^Jsi^e-^-^e^e-^' (4) 

where |0) denotes the electron vacuum in the leads, ct(i) , Ck L [t) , c pR [t) , c PR {t) are the cre- 
ation and destruction fermion operators in the leads evaluated in the Heisenberg represen- 
tation and we have made use of the identity: 

In this equation J m denotes the Bessel function of integer order m. 

We use an interaction picture to switch on adiabatically the required perturbations that 
allow us to recover the total hamiltonian H(t). In order to apply the Generalized Transfer 
Hamiltonian method we choose a hamiltonian that in its first quantized form i s§: 

H(t) = H L (t) + V L (t) = H R (t) + V R (t) (6) 
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Where Vl(£) and V R (t) can be written as: 



V L (t) = V L e vt - V AC {Q(x - x 2 ) - 0(x - x 3 )}cosu te Tlt - 2V AC Q(x - x 3 )cosu te r ' t 
V R (t) = V R e* + V AC {Q(x - x 2 ) - Q(x - x 3 )}cosu te^ + 2V AC Q(x 2 - x)cosu te^ (7) 

Vi(t) and V R {t) will be considered of the same order in the perturbative procedure and 
Vl, Vr are represented in fig 2a. The time evolution of the wave function for the total sys- 
tem can be written as: 



m)) = m E ^ 



Vac 

hu>n 



^ e -imw t e -iui kL 



Vac 

' huJa 



PR 



(8) 



This wave function must describe a particle initially on the left side. This is satisfied by 
taking /(— oo) = 1 and U R (—oo, — oo) = 0. The electrons in a particular state \ki) can in 
principle evolve to any state \p R ) in the right side so that a summation over right states 
is required in the expression of the wave function. The time evolution operator Unit, — oo) 
gives the evolution of an electron to a right state and is determined at every order from 
the time dependent Schrodinger equation by an expansion in a perturbation series. We take 
f(t) = /(°) and U R (t,—oo) = Y^jLi U R \t, — oo) where j denotes the perturbation order. 
Applying the Schrodinger equation we obtain just to order j the set of equations: 



%h 



oo BTT^ 
V J n e^t^_ e - 

dt 



PR 



\P 



R 



E Jrr 



-imujQt —iujfcj t 



V L (t)\k L 



(j) 



ihY, Jne in " ot -^e-^* t \ PR )= ]T Jn 1 e- iniUot ^e-^ t Ul- 1) V R (t)\p Rl )(j>2). (9) 

n=— oo m=— oo Pr x 

These iterative equations for Ujp are solved by projecting them on the state (p R \ ,giving : 



U R = ~ fri E E JnJm( PR \V L \k L ) 



n=— oo m=—co 
oo oo 



~ E E \[Jn+l + Jn-l]JmVAc(p R \®{x-X 2 )-<d{x-X 3 )\k L ) 

n=— oo m=— oo 



oo oo 



~ E E [ J n+1 + Jn~l]JrnV A c{PR\®(x - X 3 )\k L )} 
n=—oo m=— oo 

X /* dt ie <Vp R -Vk L -ru<>o-mwo-iri)ti 
j — oo 

For the physical parameters considered experimentally the amplitude of the time dependent 
modulation is much smaller than the barrier height, Vac << V) (Vq is the height of the 
barriers) and then it is a good approximation to consider: 
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oo oo „ t 

U R - ~* E E JnJ m {p R \V L \k L ) / dt ie *(«PH-^-™-o-™«o-^)*i (11) 

„ ^ ™_ J —oo 



n=—oo m=—oo 



Then,we proceed considering just the static part of VL(t)(V R (t)); the time evolution operator 
to j-order can be written: 



* j 1 ^ .^ u PRj-x _a,fc £ — "j-i^O— mwo-irf)tj (12) 



oo 



The solution is: 



E^ft-oo) = £££•■•• E E«» 



n,m _R X ni iij-i 



e i(u PR ~u kL -nu/ -mwo ~ijrj)t 

(u PR - u kL - nuj - mu - ijrj) 



J% 1 ~~Jn^ l (PR\VR\pR 1 }~~{pR^ 1 \VL\k L ) ^ 



(-hy (u PRi - u kL - ni^o - "T^o - «(j - l)?7)....(wp J j j ._ 1 - Wfe L - n^-i^Q - mw - irj) 
In this procedure we have assumed that the spectral density of the electrons in the leads is: 

oo 

A k {e) = 2tv JlS(e - e k + nhu ) (14) 

n=— oo 

This assumption will be correct only in the nonadiabatic regime where the external frequency 
is much larger than the inverse resonant-tunneling timecM.This regime does not consider 
the ac components of the spectral densities and Green's functions (see appendix). 
This allows us to write the Fourier transform of the retarded single electron Green's function 
of H R as: 

G + R (e)= £ J2 J n ^ ){P f , (15) 

It can be found from: 

A fc (e) = -2JmG+(e) (16) 
And then, Ujp becomes in terms of the retarded Green's function: 



oo oo p j-(ep R -E kj -nTiLu Q -mhuj —ijrj)t 
in ~ 



U$>(t,-oo) = e J2 J n J, 



n=—oo m=— oo 



(c PR ~ £k L ~ nhuj - mhuj - ijrj) 
(p R \V R G R ~(e kL + mhuo + i{j - l)r])V R ....V R G^(e kL + mhu )V L \k L ) (17) 

The transition probability from left to right per unit time can be expressed as: 
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(18) 



Where U R (t, -oo) = £°Li Ujp(t, -oo) includes the sum over all orders in perturbation the- 
ory. In the case of a continuous spectrum there are no divergences in the analytic expression 
of the Green's function so that the limit i] — > of the matrix elements in U R can be taken 
independently from that of the fractions ,and then, we can express the time evolution in 
terms of the total Green's function of the system: 

G + (e) = G+(e) + G+(e)V R G+(e) + ... (19) 

In the case of a discrete spectrum some generalization is needed because the adiabatic limit 
77 — ^ would imply to have poles in the Green's functions which coincide with these of the 
fractions that give the delta functions^. 

The expression for the time dependent transition probability is: 

p _ J T T T i(qu>o+puo-nu>o-mtjJo)t 

r RL i / j J n'J "p 1 
nmpq 

( p R \V L + V R G + (e kL +mhu )V L \k L )(k L \V L G + (e kL +qhu )V R + V L \p R y 
[ PP z t iir8(e PR - e kh - phu - qhu ) - 

PP t z iK5(e PR - e kL - nhu - mhu )] (20) 

e p R ~ e k L ~ n?luJ o ~ mhuj 

Where PP denotes the principal part and G + is the total Green's function of the system. 
In order to obtain the stationary transition probability we consider n=p and m=q: 



rl- 



n,m= 



I (PR\V L + V R G + (e kL + mhu )V L \k L }\ 2 (21) 

This formula for the transition probability is a natural extension of the Fermi's Golden Rule 
formulastirQ to analyze the coherent resonant tunneling in terms of the total Green's func- 
tion of the system. The term which does not contain the Green's function corresponds to 
first order perturbation theory, and is the only one appearing in the Transfer Hamiltonian 
methodEl, however the term containing the total Green's function is the one which includes 
processes which involve real intermediate states and therefore describes correctly the coher- 
ent resonant tunneling. 

From (21) one can obtain straightforward the Transmission probability and the coherent 
tunneling current density! . 

In many physical cases G + can be calculated but, in general, the main physics of the problem 
can be recovered by using approximations to G + . This is the case of resonant tunneling in 
which the main physics is connected with well states so that a good approximation is to 
substitute G + by the Green's function G^ of an isolated quantum well (see fig 2b) with a 
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selfenergy which accounts for the coupling of the well with the continuum of states in the 
right and left leads. By means of a Dyson equation the selfenergy can be expressed as: 

S(e fci + mhuo) = V C + V c G%(e kh + mhuj Q )V c (22) 

The imaginary part of the selfenergy is related with the elastic coupling to the leads as: 

r(e^ L + mhu ) = — 2/mE(efc L + mhu ) (23) 

And the energy shift of the resonances: 

ReJ:(e kL + mhu ) = PP f° ^ ^ = (24) 

J-oo 2-7T €k L — e qL + mnujo 

is the Hilbert Transform of the elastic coupling. 

With this procedure the energy dependence of the broadening of the resonant levels for each 
bias voltage applied across the heterostructure (this dependence is included in the potentials 
Vc) is considered in a straightforward way. 

Also its dependence on the external field is included through the dependence of £ on the 
energy of the m photo sideband (ie., the broadening of each photo sideband is different ) and 
through the spectral densities (see formula 14) that give the correct weight to each photo 
sideband. One should remark that usually the models based in the Transfer Hamiltonian 
define the elastic coupling to the leads as : 

T(e) = T L (e) + T R (e) 

= EK A ^)+J:^ R A PR (e) (25) 

k L PR 

Here T kh , T PR are the matrix elements given by BardeenEl and A kh , A PR are the spectral 
densities of the left and right leads respectively . This expression implies that the transmis- 
sion probability is the sum of the transmission probabilities of each barrier separately and 
therefore that the tunnel is sequentia We consider coherent tunnel and this is the 

reason to define the broadening as in Eq(22-23). Therefore, by means of our model we do 
not have to treat the coupling with the leads as a constant as the usual Transfer Hamiltonian 
methods for sequential tunnel do. 



III. LASER IRRADIATING THE WHOLE HETEROSTRUCTURE 

We consider now the case where the light illuminates homogeneously the whole het- 
erostructure. In this case, there is no a position-dependent phase shift of the external field 
on the heterostructure and the spatial dependence, as we will see below , appears through 
the matrix elements of the electronic momentum operator. 

The electromagnetic field is represented by a plane electromagnetic wave of wave vector k, 
parallel to the x direction and polarized in the tunnel direction (see fig.3) E = (0,0, F). 
The hamiltonian for this configuration has been solved within the framework of the Transfer 
Matrix formalism and time dependent perturbation theory.0 

Our aim now is to obtain the transition probability and the tunneling current within the 
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scheme of the Transfer Hamiltonian, in the same way as it was done in the previous section 
for the AC case. The interest of doing that in this case is again to work in terms of local- 
ized states as basis and therefore it would allow to consider, as it was discussed previously, 
different kind of excitations in the different spatial regions, as well as systems with localized 
eigenstates. 

Applying the same procedure as in the case of the oscillating voltages we can define the 
hamiltonians Hi(t) and Hn(t) including the time dependent part of the hamiltonian that 
includes only diagonal coupling between electronic states. In the Coulomb Gauge : 

H L {t) = Ek + ^ L \P z \k L ){-^—)Hae~^ + a + e^)}c + kL c kL 

Til Zt V UJq 



H R (t) = + ^(PR\P Z \PR)i^7-)H^ ot + * + e luJot )]ct R c m (26) 



jlR m* 2eVu 

Here the vector potential operator is: 

A = (—^—^(ae-™ * + a+e^'K (27) 

a + ,a are the creation and destruction operators for photons and the wave vector of the 
electromagnetic field has been neglected. We will see below that the comparison with the 
previous configuration indicates that there are qualitative features in the transmission coeffi- 
cient and current density which are very different for both cases. Therefore the configuration 
of the time dependent potential matters for analyzing the photoassisted transport proper- 
ties. 

The main reasons for those differences are, firstable, that in this case, the electron-photon 
coupling term contains the P operator matrix elements , and secondly, that there are off- 
diagonal terms in the Hamiltonian, i.e., terms which describe how the electromagnetic field 
couples the different electronic states. These terms are those which produce a difference in 
the transition probability with respect to the case without a time dependent field applied to 
the sample. In the previous case, however, the transition probability difference comes from 
the photo sidebands which appear in the regions affected by the time dependent field and 
which behave as additional tunneling channels. 

The left and right hamiltonians are exactly soluble, by means of a canonical 
transformation0y which decouples electrons and photons. We will not repit here this 
procedure but we give the expression for the retarded Green's functions associated with the 
hamiltonians Hi(t) and Hn(t): 

■ oo oo 

n=—oom=—oo 
• oo oo 

G + R (t,t') = -- £ £ Jr(ftJJ.(A a )e-^ w '^ l e- w (28) 



r=— oo s=— oo 



where: 



ft = ^k\P z \k) 
m*hujQ 
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Here, the argument of the Bessel function depends on the matrix element of the momentum 
operator and its value is very small. It is also a significative difference with respect to the 
previous case where the argument of the Bessel functions can be arbitrarily big (depending 
of the amplitude of the external field) and therefore the contribution of the terms which 
correspond to Bessel functions of higher order than zero and one, have to be considered. 
In the present case we are discussing, as the argument of the Bessel functions is close to zero 
for typical intensities and frequencies of the field , without lose of generality it is usually a 
very good approximation to neglect terms of order higher than zero or one at most. 
In order to obtain the transition probability for coherent tunneling through the double 
barrier structure, we proceed evaluating, as in the previous case , and in the framework of 
the Generalized Transfer HamiltonianEJ, the time evolution operator with: 

V L (t) = V L e* 

V R (t) = Vne* (30) 

Its expression, obtained up to infinite order in time dependent perturbation theory for the 
barriers and exact for the diagonal terms in the electron-photon coupling, is the same as 
eq.(17) but now the argument of the Bessel functions depends on the electronic momentum 
matrix elements as was already discussed. 

At this point we include the non-diagonal coupling terms in order to obtain the total time 
evolution operator. Its expression in interaction representation is : 

U(t, -oo) = 1 - l f dsH 2 (s) 
n J-oo 

+ ("i-) 2 f f drH 2 (s)U(s,T)H 2 (r) (31) 

tl J — oo J — oo 

Where H 2 is : 

eF , 

H 2(t) =J2z2^^( k i\ P *\ k i) c t' Ck i cos ( u; ot) 

^% k. 

% 

i = L,R (32) 

The caret on the H 2 operators indicates that they should be evaluated in interaction repre- 
sentation (here the "non interacting" hamiltonian is just this one we have solved exactly by 
means of the GTH) . 

We consider only first order and we keep only the Jo Bessel function terms because they are 
the only ones that give non negligible contributions due to the smallness of their arguments. 
The new term appearing in the time evolution operator has the expression: 

_ 1 J (P PR )J (p kL ) eF e ^ PR -^o-^ 

U 2 [Zj — OO J — — r 

n 2 m*ujQ (lo Pr — u kL + lo$ — irj) 

{(p R \V R Gl(e PR )P z \k L ) + (p R \P z G+(e kL )V L \k L ) + 
(p R \V R G~(e PR )V L G L (e PR )P z \k L ) + ( PR \P z G + R (e kh )V R G + (t kL )V L \k L ) 
(p R \V R G-(e PR )P z G + (e kL )V L \k L )} + {u - -u ) (33) 
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Where and G^ are the Green's functions of Hr and Hl respectively, G ± includes the 
tunneling terms and we have used the property that the total Green's function of the system 
can be expanded in terms of the two basis left and right : 

G = G L + G L V L G L + .... = G R + G R V R G R + ... (34) 

The Green's functions appearing in this expression are: 

Pr Pr ' 

GHe) = E _ ^ir— (35) 

Following the same steps as in section II the transition probability can be written as a 
Fermi's golden rule: 

Prl = -^-{\A RL \ 2 5(t PR -e kL ) + \B RL \ 2 S(e PR - e kL + hu ) + (u -> -wo)} (36) 
where A RL and contain the matrix elements: 

A RL = Ut3 PR )Ufa L )(p R \V L + V R G + (e kL )V L \k L ) 

Brl = J ° (/ H J ° (/jfcJ {(Pfll^gZ(^)^l^) + (PK|ftGS( Cjfc jVi|fc L > + 
(pR\V R G-(e PR )V L Gl(e PR )P z \k L ) + (p fl | P 2 G + R ( e^J^G+Ce^) V L |fc L > 
(pR\V R G-(e PR )P z G + (e kL )V L \k L )} (37) 

In the previous expression (37) the different terms contain the possible transitions from 
the left to the right states assisted by the light. Both resonant transitions through the well 
states and non resonant ones, for coherent tunneling are included. This kind of matrix ele- 
ments have been obtained in ref 7 in the context of sequential photoassisted tunneling. 



IV. RESULTS FOR THE TWO DIFFERENT CONFIGURATIONS 

In this section we have studied the two different configurations discussed in the previous 
section i.e., the transmission coefficient for electrons through a double barrier in the presence 
of a time dependent potential applied between the emitter and collector whose amplitude is 
position-dependent and an homogeneous electromagnetic field irradiating the whole sample . 
In fig. (4) the transmission coefficient for the first case has been drawn. The heterostructure 
consists in a double barrier of GaAs — Al x Gcii^ x As with well and barrier thicknesses of 50 
A and 100 A respectively. We have considered that the energy of the time-dependent field 
is 13.6 meV corresponding to a frequency of 3.3 THz. The important magnitude in this case 
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is the ratio between the amplitude and the energy of the field which enters in the argument 
of the Bessel functions (see formula 5). As we have already discussed depending of the value 
of this number higher order Bessel functions are non negligible and its value determines the 
number of them (i.e., the number of photo sidebands) which participate in the transmission 
probability and current density. In fig. (4) we consider that this ratio is one half and we have 
included in the calculation the Bessel functions up to fourth order. This fact is reflected in 
the four satellites which appear at both sides of the main central peak (which corresponds 
to the Bessel function of zero-order). 

One should also remark that the broadening of the resonant states in the well is not a con- 
stant but depends on the photo sideband index m (see formula 23). In this case we observe 
that the contribution of the photo sideband of index m >0 to the transmission coefficient is 
smaller than the one coming from the main peak (m=0) but of the same order of magnitude 
even for the higher order side-bands due to the finite value of the argument of the Bessel 
function Another interesting effect which was already discussed by M. Wagne is the 
quenching of the transmission probability in the presence of an AC field. This effect can 
be understood looking to formula (21) where the transition probability at the energy of the 
resonant state is weighted by : J2™=-oc Jo^n- The quenching of the transmission takes place 
if the argument of the Bessel function (3 = is such that ^0(7^) = 0. In fig (5) this 
effect is shown. We observed it for different values of (3 : as it approaches to the complete 
quenching limit the resonant transmission main peak is strongly reduced, therefore tuning 
the intensity and frequency of the time dependent modulation allows to reduce the trans- 
mission independently of the transparency of the barriers0. 

We have also plotted the current density for an AC potential through a double barrier (fig. 6) 
for different amplitudes : As main features we observe that the threshold tunneling current 
moves to lower bias and also that for higher bias the current density is smaller in the pres- 
ence of the AC field, and this effect increases as the ratio increases too. Also a step-like 
behaviour is observed in the current as a function of the external bias. We can explained 
these features in terms of the photo sidebands: the threshold, for this case is close to zero 
bias. That is due to the fact that there are photo sidebands associated to electronic states 
close to the Fermi energy which contribute to the resonant tunneling even when the resonant 
state ( E r ) is higher in energy than the emitter Fermi energy Ep. If E r is higher than Ep 
in several photon energies the photo sidebands which allow the flow of current have a low 
spectral weight and the contribution to the current is small. As far as the resonant state 
energy closes into the Fermi energy, increasing the bias, the lower indexes photo sidebands, 
i.e., those which are more intense in the spectral function can be aligned with the resonant 
state, and therefore, their contribution to the current increases. Once the resonant state 
crosses the Fermi energy the current density increases but remains smaller than the current 
with no AC applied. That is due again to the fact that the spectral function has finite 
weight in all the photo sidebands and not only in the main one which is weighted by the 
zero order Bessel function and whose value is smaller than one. In this case , only a small 
number of them , for a fixed bias tunnels resonantly and the effect of the field is to reduce 
the current. 

As the ratio between the field intensity and the frequency increases it does the coupling of 
the electron motion with the external field and the current density differs more from the 
case without this interaction and presents abrupt steps coming from the satellite bands . 
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The homogeneous illumination of the sample has been analyzed for the first time in ref 
10. In order to compare this situation with the one discussed previously, we have plotted the 
transmission coefficient for a GaAs-AlGaAs double barrier for both configurations (fig. 7). To 
make the two situations comparable we estimate the ratio supposing the total potential 
drop due to the light to be Fd where d is the length of the heterostructure, for the values 
of intensity and energy of the electromagnetic field studied we have |^ = 0.77. The main 
difference between the two configurations is that for these heterostructures the momentum 
matrix elements are very small, therefore, in the case of homogeneous light irradiating the 
sample and for typical values of the intensities and frequencies of the electromagnetic field, 
the argument of the Bessel functions of higher order than zero are negligible and therefore 
the intensity of the corresponding photo sideband can be neglected. That is the reason why 
to consider just the main peak (the m=0 term) makes sense. 

However we observe in fig. 7b the presence of two satellites in the transmission coefficient. 
These side-peaks have another origin than the photo sidebands^: They come from the mix- 
ing of electronic states due to the homogeneous light and which show up in the hamiltonian 
as the off-diagonal matrix elements of the electronic momentum coupled by the light. In this 
case just processes involving absorption and emission of one photon are considered because 
they are the main ones0. 

Therefore , the tunneling channels for the two configurations are different: in the case of an 
AC field the off- diagonal terms cancel if the amplitude of the AC potential is considered 
constant within each region (left, center and right). In this case, the main tunneling channels 
(the only ones within this approximation) are the sidebands: those in the emitter align in 
energy with the sidebands in the well producing additional contributions to the transmission 
probability and the resonant current. Their contribution can be important even for high 
order photo sidebands if the ratio is of the order or higher than one. 
In the case of an homogeneous electromagnetic field , the situation is different. The off- 
diagonal electron-photon coupling terms in the hamiltonian are those which modify the 
current density. Those channels, involving different electronic states (fig. 7 ) contribute in 
principle also with all their photo sidebands , however, as the argument of the Bessel func- 
tions ( i.e., the intensity of the m photo sideband ) is controlled by the momentum matrix 
elements and remains very small, just the zero index photo sideband (the main one) is 
non- negligible and gives a contribution to the current. Therefore the three peaks in the 
transmission coefficient come from the main bands (index zero) corresponding to three elec- 
tronic states which differ in one photon energy and which are mixed by the field. 
In fig. 8. a the current density through the double barrier in the presence of homogeneous 
light is represented. For the field intensity and frequency considered: F = 4.10 5 ^,frw7o = 
13. GmeV (3.3 THz) the effect of the external field is very small due to the electron-photon 
coupling term involving the momentum matrix elements, and the effect of the light cannot 
be observed in the characteristic curve but in the current density difference with respect to 
the case of no light present in the sample (fig. 8. b)B 

This result is in very good agreement with the experimental results of ref 3 where a double 
barrier is irradiated with a laser in the far infrared regimeP0 . 
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V. CONCLUSIONS 



In this work we have extended the Generalized Transfer Hamiltonian which allows to 
describe the coherent tunneling through a resonant heterostructurelij, to include the effect 
of a time dependent field. Two different configurations have been addressed: an homoge- 
neous electromagnetic field applied to the heterostructure and a position-dependent AC field 
dropping through the heterostructure. 

We have obtained the transition probability and the current density through a double barrier 
structure for both configurations, and the results have been discussed in terms of the main 
available tunneling channels for each situation. 

In the first case the photoassisted tunneling is mainly due to the coupling between different 
electronic states due to photon absorption and emission processes. In this case, the satellite 
bands with index higher than zero have negligible contribution to the transmission proba- 
bility. 

In the second one the sidebands at the emitter region carry the electronic charge by aligning 
with those corresponding to well and collector electronic states. The intensity and frequency 
of the AC field determines the number of effective channels (sidebands) participating in the 
current. 

We did not include the effect of the electron-electron interaction in this formalism. This 
effect is small in double barrier heterostructures but it has to be included for systems where 
the correlation is an important contribution. This is the purpose of a future work. 



VI. APPENDIX 

In this appendix we want to derive the expression for the spectral densities and to discuss 
in which limit is correct to apply the Tien and Gordon's expression for the electronic density 
of statesS. We start with the retarded Green's function which depends explicitly on two times 
t and t' but not on their difference because of the breakdown of time translational invariance 
due to the external field : 

G+(t,t') = ~e(t-t')(0\c k {t)c+(t')\0) = - l -e{t-if)exp[- % - f dtMh)] 



i g - f e k (t-f) e - 1 sinuot e j- SfcZ sinu t' 

h 



~i £ E U^U^e-^-^e-^ * (38) 

fl „,_ „ iILUq iILOq 



n=— oo m=— oo 



Let us now calculate the spectral density from the retarded Green's function which in Wigner 
coordinates r = t — t',T = is defined aslll : 

/oo T T 

Jre WT G k {T + -,T--) (39) 



and: 



G k (T + I, T - I) = l [Gi(T + I, T - I) - G;(T + I,T - I)] (40) 
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Performing the Fourier transform in the relative time we obtain: 

A k (co,T) = ^- Jn{^^)J m {^-)cos[{n - m)u T]5(uj -u k + U + m ^ ) (41) 

h n=-oom=-oc nUJ hUJ Z 

Only in the nonadiabatic regime when the frequency of the external field is larger than the 
inverse of the tunneling time for the electrons we can define an averaged spectral density as: 

1 /"To Ott 00 \l 

A k (u) = - / dTA k (co,T) = — £ J^(^)S(u -u k + ncu ) (42) 
To Jo h n= _ oc hu 

where the average is taken over one period T of the external field. 

With this procedure we restrict ourselves to the dc component of the spectral density (and 
then of the Green's functions) since the ac components (equation 41) are suppressed in the 
average. 
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FIGURES 

FIG. 1. Schematic drawing of a resonant-tunneling structure with AC voltages applied in the 
leads. 

FIG. 2. Sketch of the different Hamiltonians used to study tunneling with the GTH 
method, (a) Left and right Hamiltonians, (b) center Hamiltonian. 

FIG. 3. Schematic drawing of a resonant-tunneling structure in the presence of an electromag- 
netic field polarized in the tunnel direction. 

FIG. 4. Logio of coherent transmission coefficient as a function of energy with and without an 
AC field of parameters tyg = 0.5, hu = U.QmeV for a GaAs/AlGaAs double barrier 100-50-100 

A. 

FIG. 5. Quenching of the transmission for ^ = 2.0, = 2.4, ^ = 2.5, = 13.6meV 
(sample fig4). 

FIG. 6. Coherent tunneling current density as a function of voltage for the sample of fig4. 
(a) With and without an AC field of parameters = 0.5, Tiloo = 13.6meV. (b)For different ratios 
between the intensity and the energy of the AC field ( hujQ = 13.6meV). 

FIG. 7. Comparison of Logio of coherent transmission coefficient as a function of energy for an 
AC field and an electromagnetic field, = 0.77.(6)F = 4.10 5 ^, hio = !3.6meV 

FIG. 8. (a)Current Density in the presence of an electromagnetic field for the same parameters 
as fig7.b (b) Current density difference. 
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